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Mechanism for the high Neel temperature in SrTcC>3 



O 

(N 

O 

Q 

(N 
(N 



CO 

c3 



i 

c 

o 
o 



> 

00 

in 
in 

(N 



S. Middey 1 , Ashis Kumar Nandy 1,2 , Priya Mahadevan 2 and D.D. Sarma 3 
1 Centre for Advanced Materials, 
Indian Association for the Cultivation of Science, 

Jadavpur, Kolkata-700032, India 
2 S.N. Bose National Centre for Basic Sciences, 
ID -Block, Sector III, Salt Lake, Kolkata-700098, India. 
3 Solid State and Structural Chemistry Unit, 
Indian Institute of Science, Bangalore-560012, India 



The microscopic origin of the high Neel temperature (Tjv) observed experimentally in SrTc03 
has been examined using a combination of ab-initio electronic structure calculations and mean-field 
solutions of a multiband Hubbard model. The G-type antiferromagnetic state is found to be robust 
for a large region of parameter space, with large stabilization energies found, surprisingly, for small 
values of intraatomic exchange interaction strength as well as large bandwidths. The microscopic 
origin of this is traced to specific aspects associated with the d 3 configuration at the transition-metal 
site. Considering values of interaction strengths appropriate for SrTcOs and the corresponding 3d 
oxide SrMnOa, we find a ratio of 4:1 for the Tjv as well as magnitudes consistent with experiment. 



The picture of magnetism that has prevailed over the 
years has centred around the existence of localized elec- 
trons and their ordering leading to different types of 
magnetic order. Consequently one associates the high- 
est magnetic ordering temperatures with the more cor- 
related 3c? transition metal oxides, with examples among 
the 4c? and 5c? oxides which have wider bands, more a 
rarity than the norm. It was therefore a surprise when 
recently a high magnetic ordering temperature (T/v) of 
1023 K and 800 K were found in Ad transition metal 
oxides SrTcC>3 [l| and CaTcOs [2j respectively. These 
Tn were much higher than any of their 3d counterparts 
(SrMn0 3 ; TV=233 K) [|. In a short span, these un- 
expected experimental observations have generated a lot 
of theoretical interest and several possible explanations 
were offered [H, 0, [H . One reason was the smaller Hund's 
coupling strength and the larger bandwidth associated 
with the Ad oxides [l[. However, within a picture of itin- 
erant magnetism that one has so far, both these corre- 
spond to effects which should result in a reduction in 
the magnetic ordering temperatures and so the puzzle 
remains. An alternate explanation offered by Georges 
and coworkers was that SrTcC"3 sits at the boundary be- 
tween the itinerant to localized regime and hence has 
such a high transition temperature [5J. 

In this work we reexamine the issue of the high Neel 
temperature of SrTcOs. An unusual aspect associated 
with the half-filling of the ti g levels results in a bandgap 
opening up for a small value of U for the antiferromag- 
netic state, associated with the nesting of the fermi sur- 
face. As a result, a magnetic moment is stabilized for 
small U in the G-AFM state, while other magnetic so- 
lutions are able to sustain a magnetic moment only at 
larger values of U. In the insulating state, there are 
channels present for the electrons to delocalize and lower 



their energy only in the antiferromagnetic state. This 
energy lowering which strongly stabilizes the G-AFM 
state is larger for small intra-atomic exchange interac- 
tion strength ( Jh) as well as large hopping strength. The 
metal-insulator transition for a ferromagnetic state takes 
place at a larger value of U than the G-AFM state. So, 
although there are channels for derealization present in 
the ferromagnetic metallic state, these are smaller than 
the G-AFM state and vanish as U is increased and the 
system becomes insulating. As the energy gain by der- 
ealization for the G-AFM state decreases from 5c? to 4c? 
to 3d transition metal compounds, one expects a similar 
trend in the TV. The large Tjv observed for SrTcC-3 and 
CaTcOs we show is a generic feature of all 4c? as well as 
5c? oxides with a formal d 3 configuration on the transition 
metal atoms. 

The electronic and magnetic structure of SrTcOs has 
been calculated within a plane wave pseudopotential im- 
plementation of density functional theory using PAW po- 
tentials [|| as implemented in VASP [7]. In addition to 
the GGA form for the exchange correlation functional, we 
also included an effective U of 2 and 3 eV on Tc in the 
Dudarev implementation [8] of the GGA + U scheme as 
earlier work had shown that a value of 2.5 eV was appro- 
priate for Ru [§]. The total energies were calculated for 
different magnetic configurations using a k-mesh of 6x6x6 
k-points and a cut off energy of 400 eV for the plane wave 
basis states. In these calculations, the lattice constants 
were kept fixed at the experimental values [H, while the 
internal coordinates were optimized to minimize the total 
energy. In order to understand the origin of the observed 
magnetic stability and its dependence on microscopic pa- 
rameters, we carried out additional analysis in terms of 
a multiband Hubbard-likc Hamiltonian. The parameters 
entering the tight binding part of the Hamiltonian are 
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determined by fitting the ab-initio band structure for the 
nonmagnetic case to a tight binding model that included 
s and p states on oxygen and d states on Tc. Hopping 
is included between Tc d and O s as well as p states, as 
well as between the p states on oxygen atoms and these 
are parameterized in terms of the Slater-Roster parame- 
ters • The semi-core O s states were included to sim- 
ulate the splitting between the Tc states with t 2g and e g 
symmetry at T point H| . The Coulomb matrix elements 
entering the multiband Hubbard-like part of the Hamil- 
tonian are parametrised in terms of the Slater-Condon 
integrals F°, F 2 and F i [l2| , F 2 and F 4 were scaled 
from their atomic values to result in a particular Hund's 
coupling strength Jh , while the value of F° was chosen to 
result in an effective Coulomb interaction strength equal 
to U. A value of Jh = 0.1 eV has been used for most of 
the calculations in order to examine the small 4 limit. 
The Hamiltonian was solved using a mean-field decou- 
pling scheme for the four fermion terms till a convergence 
of 10~ 5 was achieved on the energy 13]. The total en- 
ergies determined for different magnetic configurations 
were mapped onto a Heisenberg model (-h^2J%j s,.Sj) 
with first neighbor (4) as well as second neighbor (J 2 ) 
exchange interaction strengths. An effective exchange in- 
teraction strength Jo given by 6 J\+ 12 J 2 was determined 
as this is directly related to T/v in a mean field model, 
upto a multiplicative constant. 

The total energies referenced to the nonmagnetic state 
for different magnetic configurations determined from our 
ab-initio calculations are given in Table I. As observed 
earlier [l[, the ferromagnetic and A-AFM calculations 
converge to a nonmagnetic solution for U = 0. Surpris- 
ingly the magnetic structures in which each Tc atom has 
more number of antiferromagnetic neighbors are the ones 
which converge to a magnetic solution. This was also 
found earlier and a possible reason was attributed to a 
need for beyond LDA effects due to the incorrect treat- 
ment of the residual exchange-correlation effects in Ad 
compounds Correcting for this with HSE function- 
als, they were able to examine all magnetic solutions and 
discuss trends in the Neel temperature for all members of 
the series ATcO^, where A is an alkali metal atom. We 
use a computationally less intensive method, GGA + U 
with a U = 2 and 3 eV on Tc and are able to converge 
to magnetic solutions for all configurations. Stabilization 
energies comparable with the HSE results are obtained. 
Apart from the robustness of the G-AFM ground state, 
we also find a larger stability for the G-AFM state over 
the FM state for U = 2 than for U = 3. While it is indeed 
true that the antiferromagnetic stability due to superex- 
change processes is expected to reduce as U is increased, 
it is surprising that both the magnitude of the stabiliza- 
tion energy as well as its variation with U are substantial 
for small values of U. 

The reasons for observed trends in first principles cal- 
culations are usually difficult to pinpoint as a result of 



TABLE I: Stabilization energy (meV/f.u.) with respect to 
the nonmagnetic state 



U(eV) 


Ferro 


A-AFM C-AFM G-AFM 





nonmag 


nonmag -46 -168 


2 


-432 


-568 -662 -745 


3 


-925 


-1029 -1100 -1164 



being dependent on several other parameters, which need 
not be the same as U is varied. In order to carry out a 
microscopic analysis to understand the origin of magnetic 
ordering we set up a multiband Hubbard-like model for 
SrTcOa with a U on Tc. The hopping matrix elements 
were parametrised in terms of the Slater Koster param- 
eters pda, pdn, sda, ppa and ppn [To] ]. A least squared 
error minimization procedure was used to estimate the 
best set of parameters entering the tight-binding part of 
the Hamiltonian that best fit the ab-initio band struc- 
14j . The bands with primarily Tc d character as 



ture 

well as the O p nonbonding states were included in the 
fitting. The parameters that were obtained were sda = 
-3.4 eV, pda =-3.3 eV, pdir = 1.55 eV, ppa = 0.6 eV and 
ppir = -0.15 eV and €d - e p — 2 eV. 

After obtaining the hopping strengths, we solve the 
multiband Hubbard Hamiltonian for several values of U 
and A. Features of the solutions observed earlier in the 
context of the ab-initio calculations are observed here 
also. While the G-type antiferromagnetic solution is 
found to be robust, the existence of the other magnetic 
solutions depends on the value of U. Examining the dif- 
ferent solutions at A=2 eV and a small value of U equal 
to 0.6 eV, we find that only the G-AFM solution exists 
while all other magnetic solutions converge to a nonmag- 
netic solution. Examining the Tc d and O p partial den- 
sity of states (Fig. 1(a)), we find that the t 2g states with 
primarily Tc d character are almost 3 eV wide. As a 
consequence of nesting associated with a bipartite lattice 
at half-filling [15j, a band gap opens up for just U = 0.6 
eV and .4=0.1 eV for the G-AFM configuration (Fig. 
1(b)). This increases to 1.1 eV when U is increased to 2 
eV (Fig. 1(c)), scaling almost linearly with U as shown 
in Fig. 1(d). Once the system goes insulating, it can 
very easily sustain a local magnetic moment as the hop- 
ping processes between the majority spin t2 g states at 
one site and the unoccupied minority spin t 2g states at 
the neighboring site enhances the stability of the antifer- 
romagnetic state. The variation of the magnetic moment 
in the G-AFM configuration at each Tc site is plotted as 
a function of U in Fig. 1(d). The magnetic moment is 
found to increase from a value of 0.6 \ib at U — 0.6 eV 
to an almost saturation value of 2.5 /ib at U — 3 eV, in 
contrast to a fully ionic value of 3 \ib- These values give 
a sense of the itinerant nature of the magnetism. For the 
ferromagnetic case or the other antiferromagnetic config- 
urations where some neighbouring spins are aligned fer- 
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FIG. 1: (color online) Partial density of states of Tc d states 
and O p states for (a) nonmagnetic solution with A = 2.0 
eV, U = 0.6 eV; (b) GAFM spin configuration with A = 2.0 
eV, U = 0.6 eV; (c) GAFM spin configuration with A = 2.0 
eV, U — 2.0 eV. The variation band gap and magnitude of 
magnetic moment of each Tc atom with U axe plotted in (d) 
for A = 2.0 eV 



romagnetically, the metal-insulator transition takes place 
at a larger value of U. This explains the robustness of 
the G-AFM solution in regions of the parameter space 
where the other magnetic solutions converge to nonmag- 
netic solutions. 

Two key parameters that control the properties of 
transition metal oxides are the charge transfer energy 
(A) and the onsite Coulomb interaction strength (U). 
To check the stability of the magnetic state of SrTc03, 
we have calculated the interatomic exchange interaction 
strengths J,'s with different combinations of A and U 
with J/j = 0.1 eV. The variations in J\ and J2 as a 
function of U for three values of A equal to 0, 2 and 
4 eV are shown in Fig. 2(a). As the G-AFM is found 
to be the lowest energy solution for all values of A, J1 
turns out to be antiferromagnetic while J2 (inset of Fig. 
2(a)) is primarily ferromagnetic. Considering the case of 
A=4, we find that the magnitude of J\ first increases, 
and then decreases as U is increased. The magnitude of 
J2 decreases over a small U variation and may then be 
approximated to 0. As J 2 involves higher order hopping 
processes, its value is finite only when the hopping in- 
teraction strength is sizeable (i.e. the small U regime). 
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FIG. 2: (color online) (a) Variation of Ji with U for A = 0, 2, 
4 eV. The inset shows the corresponding variation of J2. (b) 
Jo as function of U for J h = 0.3, 0.8 for SrTc0 3 and 0.8 eV 
for SrMn0 3 . (c) Color plot of J in U - A plane for SrTcOa 
with Jh = 0.1 eV. The dotted line indicates the position of 
maxima of Jo. 



As Jo is dominated by Ji, we discuss the variation in 
Ji that we find in order to understand the dependence 
of TV on microscopic interaction strengths. We are not 
able to extract J^'s in the region where only some of the 
configurations converge to magnetic solutions. A simple 
perturbative treatment of the energies was carried out 
considering a fully spin polarized ground state at each 
transition metal site and the allowed first excited state. 
This simple model is able to capture the basic physics of 
why Ji goes through a maximum as a function of U. As 
there are no derealization pathways for a ferromagnetic 
arrangement, the energy gain for the antiferromagnetic 
arrangement is directly related to J\. In the low U limit, 
the U is treated as the perturbation to the hopping and 
this tells us that Ji varies as U. In the large U limit, the 
hopping interaction is treated as a perturbation and this 
results in Ji varying as 1/U. It is immediately clear that 
this would imply that Ji should go through a maximum 
as a function of U and this is indeed what is observed 
in Fig. 2(a) for various values of A studied, all of which 
show similar dependencies. The largest J\ is seen for the 
smallest A as the effective hopping between the sites will 
be the largest there. In order to probe the role of Jh, we 
have calculated Jo for two values of Jh - a value of 0.3 
eV [lH, as well as a typical value of 0.8 eV for SrTcOa 
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and SrMn03 and plotted them in Fig. 2(b). Although 
the Jo variation shows the same trend for both values of 
Jh , there is a dramatic reduction in J in going from Jh of 
0.3 eV to 0.8 eV. This stems from the fact that the effec- 
tive exchange splitting approximately varies as U + 2 Jh- 
Similar trends are seen for SrMnC>3 for which we show 
the variation in Jo only for Jh = 0.8 eV. The parameters 
used for solving the multiband Hubbard Hamiltonian for 
SrMnC>3 were derived from a tight-binding fitting of the 
ab-initio band structure as described for SrTcC>3. The 
main difference in the extracted parameters for the two 
systems is in the value of pdir. This is 25% smaller in 
SrMn03 compared to SrTcC>3. If a U of 3 eV is believed 
to be appropriate, for SrMn03 and a slightly reduced 
value of 2 eV seems likely for SrTcC>3, the ratio of the 
calculated Jo's and therefore T^rs are in the ratio 4:1, 
with that for SrTc03 being higher and consistent with 
experiment. Using the caculatcd Jo for SrTcC>3 and as- 
suming a mean-field scaling of 80%, we get a Tjy of 1020 
K, in very good agreement with experiment. The values 
of Jo have been extracted in the complete U - A plane, 
and plotted in Fig. 2(c). The largest values of Jo are seen 
near U = 3 eV for small A (~ 0.3 eV) and as a function 
of U are found about the dashed line drawn with a slight 
slope to the line U = 3. There are two dominant ener- 
getics which determine the magnitude of Jo- The first is 
the delocalization channels present in the antiferromag- 
netic arrangement and the second are the delocalization 
channels present in the ferromagnetic metallic state. As 
a result of nesting, while the antiferromagnetic solution 
goes insulating at small values of U, 0.6 eV for A = 2.0 
eV, the ferromagnetic solution remains metallic upto a 
value of U equal to 2.5 eV which is of the order of the 
bandwidth. Hence depending on the magnitude of the 
effective hopping parameter, one could be in a regime 
where the peak is close to the point where the ferro- 
magnetic solution goes insulating or far from it. With 
increasing A as the effective hopping strength between 
neighboring Tc site decreases, one finds a reduction in 
Ji. Additionaly for a fixed value of A, as U is increased 
one has a maximum and then a decrease with all Jo's 
merging to the same value at large U when the hopping 
is too small to bring about any magnetic ordering. 

The results of Fig. 2 indicate a wide parameter regime 
for the stability of the G-AFM state as the ground state, 
thereby suggesting that the large magnetic stabilization 
energies are not limited to just SrTcC>3 and CaTcC>3, 
but is generic of all transition metal oxides with a d 3 
configuration. Another probe which could be used in ex- 
periments is hydrostatic pressure which could be used to 
examine the stability of the G-AFM ground state. Evalu- 
ating the total energy differences between the G- and the 
C-AFM ground states within ab-initio GGA based cal- 
culations using the SrTc03 structure as a starting point, 
we find that with a small change in the lattice constant 
of 2.5%, the C-AFM solution is no longer stable and con- 



<3 -50 

■u -100 

S -150 

< -200 



- - » 



(a) 



GAFM-Nonmag 
-CAFM-Nonmag 



5 10 15 

% decrease of each lattice constant 



\P) Expt. structure 


Tc d partial dos 

O p partial dos 


llnwn cnin 


ft p=-^ 

i i i i i i 


— 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 

■ fojtfb. ' O rrf 


Up spin 


■ 5% decrease of each lattice 
- constant 


U" Down spin 

H-H-t-H 


. i 1 I 1 i 1 i 1 I 1 i 


Up spin 

A — C>*» 


• 10% decrease of each lattice 
- constant| , ■ , | , | , 


Vr Down spin 
i.i.i. 




V 


1 1 1 1 1 1 

Up spin i 

h— W 


15% oecrease of each itf 1 
lattice, constant , i 


<\ 


Down spin | 



-12 -10 -8 



-4 -2 2 4 
Energy (eV) 



FIG. 3: (color online) (a) Variation of magnetic stabilization 
energy for GAFM and CAFM spin confiurations with the per- 
centage decrease of each lattice constant, (b) Tc d and O p 
partial density of states for GAFM configurations for different 
compressive strain. For 15% decrease of lattice constant, the 
final spin configuration is nonmagnetic. 



verges to a nonmagnetic solution (Fig. 3(a)). The G- 
AFM solution we find is stable and as revealed by the 
total energy difference with the nonmagnetic state plot- 
ted in Fig. 3(a) as well as the Tc d and O p projected 
partial density of states plotted in Fig. 3 (b) for 0, 5, 10 
and 14% decrease in the lattice parameter. For a lattice 
parameter change larger than 10 % we find a collapse of 
the G-AFM solution to the nonmagnetic solution, again 
supporting the robustness of the G-AFM solution. 

We have examined the origin of a high Tjv in SrTc03 
within the mean-field limit of a multiband Hubbard 
model. A wide range of parameters is found for which 
the G-AFM solution is stable. In contrast to usual ex- 
pectations, the magnetic stabilization energy is found to 
be larger in the limit of small Jh and larger bandwidths. 
This is traced to aspects of the d 3 configuration at the 
transition metal site, which in the insulating state, al- 
lows electrons to delocalize only in the antiferromagnetic 
configuration. Using the appropriate values of the in- 
teraction strengths for SrTc03, as well as a mean- field 
reduction factor, we get a Tn of 1020 K, which is four 
times larger than that for SrMnOs. 
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